Optimal quantum estimation of loss in bosonic channels 
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We address the estimation of the loss parameter of a bosonic channel probed by Gaussian signals. We derive 
the ultimate quantum bound on precision and show that no improvement may be obtained by having access to 
the environment degrees of freedom. We found that, for small losses, the variance of the optimal estimator is 
proportional to the loss parameter itself, a result that represents a qualitative improvement over the shot noise 
limit. An observable based on the symmetric logarithmic derivative is derived, which attains the ultimate bound 
and may be implemented using Gaussian operations and photon counting. 
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In all the branches of physics, many quantities of interest 
are not directly accessible, either in principle like the mea- 
surement of all fields [1] or due to experimental impediments. 
In these cases, one should resort to indirect measurements, 
inferring the value of the quantity of interest from its influ- 
ence on a given probe. This is basically a parameter estima- 
tion problem whose solution, i.e the determination of the most 
precise estimator, unavoidably involves an optimization pro- 
cedure [2]. Indeed, the central result in classical estimation 
is the so-called Cramer-Rao inequality, which sets a lower 
bound on the variance of any estimator in terms of the Fisher 
information. When quantum systems are involved, the op- 
timal measurement to detect an unobservable quantity may 
be found using tools from quantum estimation theory J51 01 • 
The quantum version of Cramer-Rao inequality has been es- 
tablished and the lower bound imposed by the 
quantum Fisher information has been shown to be achievable 
asymptotically |6]. In turn, this results permit to write a pa- 
rameter based uncertainty relation for unobservable quantities 

EH. 

Quantum estimation theory has been mostly applied to find 
optimal measurements and, in turn, to evaluate the corre- 
sponding lower bounds on precision, for the estimation of 
parameters imposed by unitary transformations. For bosonic 
systems these include sing le-mode phase ifiHllfll, displace- 
ment fllll . squeezing 1 12, 3 as well as two-mode transfor- 
mations, e.g. bilinear coupling 1 1-411 - Concerning open quan- 
tum systems and non unitary processes, QET has been ap- 
plied only to finite dimensional systems [15], to optimally es- 
timate the noise parameter of depolarizing 1 1611 or amplitude- 
damping IU7I1 channels. On the other hand, to the best of our 
knowledge, QET has been not exploited to determine the opti- 
mal way of estimating parameters of a lossy bosonic channel. 
Needless to say, besides fundamental interest, quantum lim- 
its to the estimation of bosonic lossy channels are extremely 
relevant for applications, as for example absorption measure- 
ments and characterization of optical media |18, 20ll - 

The precision of a measurement corresponds to the smallest 
value of the parameter that can be discriminated. In several 
quantum mechanical schemes, precision is improved by the 
more efficient use of the resources employed in the measure- 



ment process 112 ll 12211 . A relevant example is the quadratic in- 
crease of the precision in phase estimation by squeezed states. 
On the other hand, despite the improvement in the resource- 
precision trade-off, known estimation schemes are character- 
ized by a threshold value, under which the parameter cannot 
be measured. The results described in this letter, besides being 
the first example of optimized parameter estimations for non 
unitary maps, show a novel feature of quantum-limited mea- 
surements: the asymptotic variance scales with the parame- 
ter itself. Therefore, our results may be applied to arbitrarily 
small values of the parameter, with a fixed relative variance. 

The scheme we are going to consider is the following: the 
dynamics of a bosonic quantum system is governed by a Lind- 
blad Master equation of the form 
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C[a]g . 



(1) 



which results from the interaction of the system with an ex- 
ternal environment, as for example a bath of oscillators. We 
want to estimate the value of the loss parameter 7, and to this 
purpose a set of identically prepared signals (probe state) are 
sent through the channel and measured at the output. We de- 
note by go the input state and by g 7 = E^(qq) the state at 
the output of the channel, E 1 being the map associated to the 
evolution (fl}. An estimation strategy for 7 consists of three 
ingredients: a suitable input signal, a measurement at the out- 
put and an inference rule to extract the value of the parameter 
from the experimental sample x, i.e. the set of measurement 
outcomes. The ultimate goal of choosing an estimation strat- 
egy is to find a scheme that maximizes precision with the con- 
straint of a finite fixed amount of energy impinged into the 
channel. In particular, in this letter we address the following 
questions: 1) Which is the best probe state? 2) Which is the 
optimal measurement that should be performed at the output? 
3) Which is the attainable precision? 4) Can the precision be 
improved by accessing the environment degrees of freedom, 
e.g. the bath of oscillators? 

The quantum state at the output of the channel may be 
formally written as g 1 = E 1 (go) = exp (^t/2C[a\) go, 
whereas the constraint of finite probe energy reads as follows 
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t„.l = 



it] = 1 being the bosonic mode operators. 
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In this letter we focus attention on Gaussian probe states and 
consider a zero-temperature environment without squeezing, 
i.e. we focus on lossy channels described by a superoperator 
of the form 



£[ci]f? = 2a) ga — a)ag - 



ga^a . 



(2) 



Given the evolution time (length) t, the aim is to estimate as 
precisely as possible the parameter 7 in the Lindblad Master 
equation ((T). In the following we consider the parametriza- 
tion exp(— yt) — cos 2 <fi and seek for the best estimation 
strategy for tj>. We denote the probe state as go and the 
evolved state as g^ with the Master equation that rewrites as 
dg/dcf) = tan <j)£[a] g. Our goal is to determine the optimal 
probe state go and the optimal measurement, i.e the proba- 
bility operator-value measure (POVM) O x in order to form 
the estimator <fi to infer the the value of <fi from the set of 
outcomes \- Notice that, being our estimation problem uni- 
parametric, the optimal measurement can be realized using in- 
dividual measurements on separate (subsequent) preparations 
ll23ll . so that the expectation value of the estimator can be writ- 
ten as E [</>] = J2 X 4>{X) T±[e<t>)® N O x }. We focus attention 
on asymptotically unbiased estimators, i.e. those for which 
E0[0] — > 4> as N — > 00. Those are interesting because as the 
number of experimental data increases, all systematic errors 
go to zero. The variance of an unbiased estimator is bounded 
by the Cramer-Rao inequality 02], 



Var [0] > 



1 



(3) 



where F(<f>) is the Fisher information of the measurement. 
F((j)) is additive, i.e the total Fisher information of measure- 
ments on multiple copies of g^ is the sum of the individual 
ones. It is worth mentioning that the so-called maximum like- 
lihood estimator asymptotically attains this bound for large 
number of identical repeated samplings of a probability dis- 
tribution. In such case the asymptotic variance scales as 1/N. 
i.e. Var [0] -> [NF (</>)]- x . The Fisher Information F(0) 
is bounded from above by the Quantum Fisher Information 
(QFI) F(<f>) < H(4>) which, in turn, provides a measure of 
the ultimate precision available with a given quantum state. 
QFI is additive too and thus we have 



Var, 



> 



1 



NH{(f>) 



(4) 



where H ((f)) can be expressed in terms of the symmetric loga- 
rithmic derivative, (SLD) A(cf>). SLD is implicitly defined as 
the Hermitian operator that satisfies 



dgcj, 



1 



[^A(0) + A(<f>)g^] 



(5) 



from which the QFI can be computed as H((f>) = 
Tr[g0 A(<fr) 2 }. A(<fi) is an observable which depends on (f> and, 
analogously to the classical LD, its expectation value is zero, 
i.e. Tr[A(0) g^] = 0. A biased A(cj)) reveals to which extent 



the measured state differs from g^. In addition, it has been 
proved [5] that an optimal measurement can be obtained by 
projecting onto one-dimensional eigenspaces of A(0), i.e. a 
measurement for which the Fisher information is maximized, 
F(cj>) = H(<f)). Upon this considerations, in the following we 
consider the one-step adaptive strategy EB, in which one 
makes a rough estimate <fio on a vanishing fraction of copies 
N s , 1/2 < S < 1, and then measures A(</>o) on the remaining 
copies in order to refine the estimate (p. 

Notice that Eq. © can be regarded as a Heisenberg rela- 
tion for parameter estimation. In the case of pure states and 
unitary evolution of the form U(<fi) = exp(iG0) we have 
H(4>) = 4(AG) 2 , where (AG) 2 is the uncertainty of the 
generator G on the state used as a probe. Overall, we get 
the inequality Var0[0](AG) 2 > 1/(4N), which represents a 
parameter-based uncertainty relation for unobservable quanti- 
ties flU. 

Any Gaussian state of a single bosonic mode may be 
represented as a thermal state p^ under the action of 
a squeezing and a displacement operations i.e. g — 
D{a)S(C,)p^(C)D^(a), where a = se w is the displace- 
ment amplitude and £ = re~ 2%v is the squeezing parameter; 
D(a) = exp(aat - a*a), S*(C) = exp^V 2 - iC* 2 a 2 )- 
The thermal state can be expressed in terms of its purity as 
= 2fi/(l+fi)[(l~fi)/(l+fi)] a '' a . At this point, some con- 
siderations about the choice of a probe state are in order. We 
aim at finding the Gaussian state g which, under the action 
of the amplitude damping channel, is mapped onto g,p with 
highest QFI. Of course, if an arbitrary amount of energy were 
available, an infinite precision could be reached. If instead 
one restricts to a probe with finite energy, then a compromise 
between \a\, r and p, must be achieved. From the phase sym- 
metry of the amplitude damping channel, it is clear that the 
only relevant angular parameter is the difference between the 
displacement phase arg a and the squeezing direction arg £. 
Hence, we can take (et£K and consider 9 = arg a as the 
relevant angle. On the other hand, as also intuitively expected, 
it is of no use to spend energy in preparing a thermal state, 
that is, for any given input energy n, the optimal probe state 
is pure. Therefore, the optimization problem can be reduced 
to determine two parameters, the the ratio x of squeezing en- 
ergy to total energy and the displacement phase 9. Since for 
a pure Gaussian state go = D(ao)S(r o )\0) (0\S^ (r )D^ (ao) 
the mean photon number is given by n — sinh 2 ro + \ao\ 2 we 
have sinh 2 ro = xn and |ao| 2 = (1 — x)n. 

The evolution of the state parameters under the action of 
the channel may be explicitly evaluated [24] 



p^ = 1/ y cos 4 <p + sin 4 4> + 2 cos 2 <fi sin 2 <fi cosh 2ro 
= i cosh -1 [p<p (cos 2 4> cosh 2r + sin 2 ^) ] 



S4> 



so cos ( 



(6) 



We can now evaluate the SLD A(</>), which corresponds to the 
optimal measurement and allows to calculate the QFI. Upon 
writing g^ in its diagonal form, = J2k Qk\^k){^k\, \i>k) = 



3 



D(a)S(r)\k) one easily finds 

{ll>q\£[a]Q\lpp) 

pq 



A((t>) = 2tan0 



Qp + Qq 

A lengthy but straightforward calculation yields 



\lPq)(lPp\ 



A(4>) = 2tan(/)D(a)S(r)KS f (r)D^ (a) 
K = [Aa)a + B{a 2 



(7) 



(8) 
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where 

A = 



2 ft 



1-p 



-(^cosh2r-l) B 



/'" 



1 + M 2 



■ sinh 2r 



C = fi(a cosh r + a* sinh r) F = 1 



2/xcosh r 
1 + M 2 



(9) 



All the parameters in (O are given in Eq. (0, i.e. refer to the 
state after the evolution in the lossy channel. We omitted 
the explicit dependence on <fi for brevity. Using the identity 

K = F' + tS\t])D(/3) a?aD(/3)S(<n) (10) 

where F' = F + ±(r - A - 2r\/3\ 2 ), tan2f? = -2B/A, 
t = VA 2 - AB 2 and (3 = (Ccoshr + C* sinh r) //rr, 
one sees that the eigenvectors of A(</>) are of the form 
D(a)S(r)S' i [rj)D^ (f3)\n), which means that the measure- 
ment of A(0) may be implemented with Gaussian operations 
and photon counting, e.g. by squeezing and displacing the 
state under investigation and then measuring the photon num- 
ber distribution by a suitable reconstruction technique [25]. 
The explicit evaluation of the QFI H(<fi) yields 



Azn 



1 + z{2 + z + Anx) 

Anx 2 z(l + fix) 
1 + z(2 + z + 2nx) 

alt 



l- x + 2nx^ Vz (11) 

z 



+ 2(1 - x)yjnx(l +nx) 



where z = e 1 — 1 = tan 4>, and where we have already 
performed the trivial optimization over the phase 9, which 
yields 9 — 0, i.e. the displacement should be performed along 
the same direction of squeezing. The optimization procedure 
thus reduces to maximizing the QFI H (cf>) with respect to the 
squeezing ratio x. In Fig. (HJ we report the renormalized QFI 
H{(f))/fi as a function of the squeezing ratio for different val- 
ues of n and of the actual loss parameter. As it is apparent 
from the plots in the regime of small losses and small probe 
energies the optimal probe is the squeezed vacuum (x opt = 1), 
whereas for increasing energy there is a nonzero value of the 
optimal squeezing fraction x opt , which is a monotonically de- 
creasing function of both the probe energy and the loss param- 
eter itself. In the small energy regime, with squeezed vacuum 
as the optimal probe, the QFI reads 



H{4>) = T 



4n(l + z 2 ) 



+ 2z(l 



4n + 0< 



(12) 



where the second equality expresses the attainable precision 
in estimating the loss parameter for weakly damping channels. 
In the following, we will see that this is the ultimate limit even 
if one has access to the environment degree of freedom. In 
terms of 7 the bound reads 



Var 7 [ 7 ] 



7 



nNt 



Oil 2 



(13) 



This is a remarkable result since it is valid for any value of the 
loss parameter with no lower bound (see also Fig. [2j. Recall, 
however, that squeezed vacuum probes will not be optimal for 
large enough values of 7. 



H(<*)/ n 




FIG. 1: Normalized Quantum Fisher information H (<f>) /n as a func- 
tion of the squeezing fraction of the probe state for different values of 
the probe energy and two values of the actual loss parameter. (Left): 
tan 2 4> = 0.1. (Right): tan 2 (f> = 5. In both plots, from bottom to 
top in the region x ~ 0, the curves for n — 0.5, n = 1, n = 2, 
n = 5, n = 10, and n = 100. 

The maximization of the QFI in the general case may be 
done numerically. In Fig. [2] (left) we report the log-log plot 
of the rescaled optimal variance AVar^,[0] = 1 / H max (4>) as 
a function of the probe energy n for different values of the 
actual loss parameter. As it is apparent from the plot the vari- 
ance does not dramatically depend on the actual value of the 
loss parameter. The common scaling is given by Var^, [</>] oc 
(nN)' 1 . 
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FIG. 2: (Left): Rescaled optimal variance H~ 1 (cj>) as a function of 
the probe energy for different actual values of the loss parameter: 
z = 0.01, 0.1, 1, 5, 10. The larger is z (<j)) the closer is the curve 
to the asymptotic _ff _1 (<j!>) = (4n) _1 . (Right): Rescaled optimal 
variance compared to the variance obtained for coherent probe as a 
function of the loss parameter for different values of the probe en- 
ergy: from top to bottom the curves for n = 0.1, 0.3, 0.7, 1 respec- 
tively. The improvement obtained using optimized squeezed states is 
apparent. 

In order to better understand the behavior in the large en- 
ergy regime an asymptotic analysis is in order. Upon expand- 
ing the optimization equation dH/dx = around z = and 
retaining the leading order we obtain analytically the op- 
timal value of x for small 7, which for large n can expressed 



4 



as x op t = {Anz)~ x / 2 + 0{l/n). This shows that, when a 
large amount of energy is available, one can improve perfor- 
mances, besides squeezing, by employing a fraction of en- 
ergy to displace the probe. The QFI in this regime reads 
H{4>) = 4z _1 (1 — y/nz) + (2 + 4n) + 2z, which corresponds 
to a variance Var^f^] — > j- (1 _|_ ^n<j>) + O (^ 4 ) Finally, by 
rewriting this in terms of 7 we have 

^^■m +o {yfw)^ (i4) 

that is, the proportionality of the variance to the parameter is 
recovered for large n. Notice that the scaling H(<j>) oc n is 
valid also in the general case, as it is apparent from Fig. [2] 
(left). 

It is worth stressing the improvement of precision with re- 
spect to that attainable with coherent states (shot-noise limit). 
In fact, by taking x = one easily sees that in this case the 
QFI reads H((f) — 4nz/(l + z) which corresponds to to 
Var 7 [7] — > (nNt 2 ) -1 . In other words, the proportionality 
of the variance to 7 cannot be achieved using coherent probe 
states. In Fig|2](right) we report the optimal rescaled variance 
as a function of the loss parameter compared to the 
variance that can be obtained using coherent probes. The im- 
provement at small values of the loss parameter is apparent. 
This can be intuitively understood as follows: for small val- 
ues of 7 the action of the loss map E 1 on a coherent probe is 
only that of a small displacement, whereas a squeezed vacuum 
is dramatically "mixed up". Hence, a squeezed state is more 
sensitive to small losses. Contrarily, for large values of 7 a 
coherent state is highly displaced whereas a squeezed vacuum 
state becomes close to pure again, losing its sensitivity to 7. 

Let us now discuss other ways in which this performance 
could, in principle, be improved. The Heisenberg limit gen- 
erally describes the ultimate precision attainable for parame- 
ter estimation. As we have seen above, this holds in a strict 
sense when one deals with unitary transformations, where 
the QFI characterizes the sensitivity of a state for the esti- 
mation of a parameter. However, when one deals with non- 
unitary transformations, this may not be true because the 
probe state evolves into an entangled state with the environ- 
ment. The access to the degree of freedom of the environ- 
ment may provide improved precision to the measurement. In 
the case of the loss of a channel, some of the energy present 
in the probe state is lost through coupling to the environ- 
ment. The Master Equation (1) can be seen as the effec- 
tive interaction of the mode a with a second vacuum mode 
b through a bilinear (beam-splitter-like) evolution of the form 
U{<j>) = cxp \i<j) (a)b + ab^)] . This allows to have a uni- 
tary representation of the process. Here we use this picture 
to derive the attainable precision in the hypothetic case that 
one had access to the bath of oscillators. In such a situa- 
tion, the state under inspection would remain pure, therefore 
A(^) = 2dp/d(f>, as can be seen by taking the derivative of the 
identity p 2 = p. The generator G is given by G = (a^b+atf), 
and the uncertainty (AG) 2 is (AG) 2 = n. As a consequence 



one gets Var^ \(f>] — > , which corresponds to the precision 
attained using squeezed vacuum probe. 

In conclusion, we have shown that using Gaussian squeezed 
probes one can improve the estimation of the loss parameter 
of a bosonic noisy channel. As it holds for any single parame- 
ter quantum estimation problems, a one- step adaptive scheme 
is optimal to leading order, thus providing a practical way for 
implementation. We have shown that the optimal measure- 
ment for Gaussian probe states can be implemented by means 
of Gaussian operations and photon counting. Furthermore, 
for small losses, the estimator variance obtained by choosing 
appropriate probe states decreases proportionally to the loss 
parameter, thus providing unlimited resolution for arbitrary 
small losses. We have also obtained the optimal trade-off be- 
tween squeezing and displacement of probe states, showing 
that squeezed vacuum states are optimal in the small energy 
limit. Finally, we have shown that even by having access to 
the environment, one cannot improve the performance of our 
scheme. 
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